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Abstract 

We present a preliminary determination of the CKM matrix element \Vub\ based on the analysis of 
semileptonic B decays from a sample of 88 million T(4S') decays collected with the BABAR detector 
at the PEP-II e~^e~ storage ring. Charmless semileptonic B decays are selected using the electron 
energy E^. and the invariant mass q'^ of the electron-neutrino pair. The neutrino momentum is 
inferred from a measurement of the visible energy and momentum in the detector and knowledge 
of the e^e~ beam momenta. The partial branching fraction is determined in a region of the 
q^-Ef, plane where semileptonic B decays to charm are highly suppressed. The total charmless 
semileptonic branching fraction is extracted using a theoretical calculation based on the heavy 
quark expansion. Preliminary results yield |Kfe| = (4.57 ± 0.21 ± 0.25 ± 0.34 IgJ^ ± 0.22) x lO^^ 
where the uncertainties are from statistics (data and MC), detector modeling, background modeling, 
the shape function, and the heavy quark operator product expansion, respectively. 
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1 INTRODUCTION 



The study of the weak interactions of quarks has played a crucial role in the development of the 
Standard Model (SM), which embodies our understanding of the fundamental interactions. The 
increasingly precise measurements of CP-violating asymmetries in B decays allow stringent exper- 
imental tests of the SM mechanism for CP violation [1] via the non-trivial phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. Improved determinations of \Vub\, the coupling strength of 
the b quark to the u quark, will improve the sensitivity of experimental tests of the SM description 
of CP violation. 

Early determinations [2] of \Vub\ were based on the study of the lepton momentum spectrum 
near the kinematic endpoint in semileptonic B decays, where the background from the dominant 
decay chain® b ceV is kinematically forbidden. See Ref. [3] for more recent measurements using 
this method. This technique provides precise measurements of the partial branching fraction for 
lepton momenta^ above ~ 2.3 GeV/c, but accepts only ~ 10% of the b — ^ ueV rate, resulting in a 
substantial theoretical uncertainty in the determination of \ Vub\- The major part of this uncertainty 
comes from limited knowledge of the "shape function" (SF), i.e. the distribution of the b quark 
momentum inside the B meson [4]. These uncertainties can be reduced by measuring the energy 
spectrum of photons from the decay 6 — > 57 [4, 5], but this is very challenging [6]. 

Recently, measurements of \Vub\ that use the invariant mass mx of the hadronic system in 
semileptonic B decays have appeared [7, 8]. The use of mx to select b ueu decays results in a 
much higher acceptance (~ 70%) of the decay rate. It is, however, experimentally challenging, since 
the association of particles with the semileptonic B decay is rendered difficult by the presence of the 
decay products of the B in the event. These measurements also have significant SF uncertainties. 

In this analysis a new approach [9] is taken to the determination of \Vub\- Semileptonic B 
decays are selected using energetic electrons and simultaneously making requirements on q^, the 
invariant mass squared of the eV pair. The neutrino 4-momentum is reconstructed from the visible 
4- momentum and knowledge of the e"'"e~ initial state. The dominant charm background is then 
suppressed by selecting a region of the — E(, phase space where properly reconstructed b ceV 
events are kinematically excluded. The fraction of the b ueU phase space accepted is ~ 20%. The 
amount of background events remaining in the signal region due to resolution effects is evaluated 
in Monte Carlo simulations. The determination of \Vub\ in this method is sensitive to the b quark 
mass through both the heavy quark operator product expansion [10] (HQE) and the SF, but has 
little sensitivity to Fermi motion as described below. 

2 DATASET AND SIMULATION 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy e+e^ storage ring between 2000 and 2002. The BABAR detector is described in detail 
elsewhere [11]. A sample of 81.4 fb~^ collected at the T{4S) resonance, corresponding to 88.4 
million BB pairs, is used along with 9.6 fb~^ collected at center-of-mass energies approximately 
40 McV below BB threshold. The data below BB threshold, scaled in cross-section and luminosity, 
and whose particles are scaled in energy, are used to subtract the non-BB contributions from the 
data collected on the T(45) resonance. Simulated BB events are used in estimating efficiencies 
and backgrounds. Branching fractions and form factors are taken from Ref. [12] in most cases. 

^Throughout this paper, whenever a mode is given, the charge conjugate is also imphed. 
'^Throughout this paper, quantities are given in the T(4S') rest frame unless stated otherwise. 
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Branching fractions for the semileptonic B decays to charm are adjusted as described below. The 
simulation of charmless semileptonic B decays is based on the model described in Ref. [13], which 
calculates the triple differential decay rate to order as and convolutes it with a shape function 
parameterized as 

F(A;+) = Ar(l-xf e(^+")^; x= < 1, 

TUB — 

where k+ is the residual 6-quark momentum, to account for the non-perturbative interactions of 
the b quark within the B meson. We use [14] mf^ = 4.735 GeV/c^ and a = 1.6. This model 
produces a spectrum of hadronic masses mx down to 2mj^ that does not contain any resonant 
states. Subsequent fragmentation of the mesons was simulated via JETSET [15]. Decays to low- 
mass hadrons (tt, r), p, uj, r)') are simulated separately, using the form factor model of Ref. [16], 
and mixed with the non-resonant states in such a way as to keep the mx, and Eg spectral 
distributions the same as in the inclusive model. 



3 ANALYSIS METHOD 

Events are selected based on the presence of an identified electron with energy Eg > 2 GeV using 
the criteria given in Ref. [17]. Electrons from the decay J/iJj e+e~ are vetoed. The criteria 
given in Ref. [17] for selecting hadronic events and rejecting radiative Bhabha events are applied. 
Within these events the total visible 4-momentum pvis is determined using charged tracks emanating 
from the collision point, identified pairs of charged tracks from Kg — > tt+tt", A — > p7T~ and 
7 e^e~, and energy deposits in the electromagnetic calorimeter. Each charged particle is 
assigned a mass hypothesis based on particle identification information. Only those calorimeter 
clusters unassociated with a charged track and which have a lateral energy spread consistent with 
photons are considered. Energy deposits due to long-lived neutral hadrons are not efficiently 
rejected by these criteria and are treated as coming from massless particles in the 4-vector sum. 
We form the missing 4-momentum Pmiss = Pe+e- ~ Pvis, where Pe+e- is the 4-vector of the initial 
state. The 4-vector pi, = (Pi,,\P^\) is used as an estimate of the momentum of the neutrino 
from the decay b ucu J wlicrc is derived, froni Pmiss 

by applying a bias correction that was 
determined from Monte Carlo signal events. Additional requirements are made to improve the 
quality of the neutrino reconstruction and suppress contributions from e"'"e~ qq continuum 
events. Each event must satisfy (1) no additional identified leptons, (2) —0.95 < cos^miss < 0-8, 
(3) 0.0 < ^miss - jPrnissjc < 0.8 GcV, (4) 0.0 < jPmissI < 2.5 GeV/c and (5) |Pe • T| < 0.75 |Pe| |T|, 
where ^miss is the angle of the missing momentum vector with respect to the beam axis, and T is 
the thrust vector for the event, excluding the signal electron candidate. The large background from 
b ceV decays is then suppressed by calculating = {Pe + Pv)"^ and computing the maximum 
kinematically allowed hadronic mass squared, s™^, for a given Eg and q^. In the case where 

±Ee > i-^'y- (t^)' maximum invariant hadronic mass squared is 




otherwise 



where (3 = 0.06 is the boost of the B meson in the T{4S) frame. We require < 3.5 GeV^/c"^ 2± 
mjj. Figure 1 shows the distribution of generated b — > ceV and b — > ueU decays in the q^-Eg plane 
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and indicates the contour corresponding to s^'"^^ = m^). The requirements for E^, and s™^^ as well 
as criteria (l)-(5) were optimized to minimize the total (experimental and theoretical) uncertainty 
on B{b ueV). 




Figure 1: Distribution of versus Ef, for generated h — > cev (left) and h — > uev (right) events 
in the T{AS) rest frame. The curved contour corresponds to s™^^ = and the diagonal line to 
= GeV^/cl 



The quality of the neutrino reconstruction was evaluated using a control sample (DeV) consist- 
ing of ~ 90, 000 decays of the type B —>■ D^eV{X) where the is reconstructed in the K~t:^ decay 
mode and satisfies \PdA > O-^ GeV/c and the electron satisfies Ee > 1.4 GeV. The D^-e combina- 
tion must satisfy —2.5 < cosOb-Db < 1-1 where cosOb-Ds = C^EbEd^ — m?B ~ "^"De) li'A^ b\\P De\) 
is the cosine of the angle between the vector momenta of the B and the D^-e system under the 
assumption that the only missing particle in the B decay is a single neutrino. This criterion selects 
both B — > D^eV and B — > D*eV, D* — > D^{7r, 7) decays with high efficiency while rejecting 
other sources of D^-e combinations. The additional requirements made on the semileptonic B 
decay in the control sample selection differ from the signal selection, where only the electron is 
required. However, in neither case are any requirements made on the decay of the other B, leaving 
its properties the same in both samples. An estimate of the neutrino energy can be formed from 
the known B energy and the measured and electron energies. A second estimate of the neutrino 
energy is given by the \Fu\ defined above. In simulation the first (second) estimate for the true 
neutrino energy has a bias of 0.014 (0.271) GeV and an r.m.s. of 0.202 (0.359) GeV for Deu events. 
Subtracting the first estimate from the second gives the distribution shown in Figure 2, from which 
we find a mean (sigma) of 0.109 (0.433) GeV on data and 0.107 (0.427) GeV on MC, indicating 
that the missing energy from the other B d.eca.y is weU modeled. The distribution of |Pmiss 

I, also 

shown in Figure 2, is sensitive to both the modeling of the semileptonic b — > ceU decays and of 
missing energy, and shows good agreement with the simulation. 

The Deu control sample was also used to improve the modeling of the b ceU decays as 
follows. A binned fit in the variables \Pd\i Eg and cos 6 b-Dc ^ is performed on the DeV sample 
after subtracting continuum and combinatorial background. The fit determines scale factors for 
the MC components B — > DeV, B D*eV and other contributions, of which 85% are B D**eV 

*The cosOb-Db requirements were relaxed to —10 < cosOb-Db < 5. 
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Figure 2: The top plot shows the difference between the two neutrino energy estimates described in 
the text for continuum-subtracted data and BB MC. The bottom plot shows the |Pmiss| distribution 
for the DeV control sample. 



(P-wave charm mesons) or i? — > D^*^7reV; the remainder comes from B DD decays and from 
misidentified electrons. Updated branching fractions are obtained by requiring the total b — > ceU 
branching fraction to equal the measured value [18]. The fit raises the branching fractions relative 
to those of Ref. [12] for B — > Del7 and B D*eV and lowers the remaining contributions. Given 
the large uncertainty in the composition of the high mass charm states in semileptonic B decay, 
this procedure cannot be considered a measurement of the branching fractions, but rather a means 
of improving the simulation of inclusive quantities in 6 ^ ceV decays. We find an improved 
agreement between data and MC in kinematic distributions {Eg,, Pmissj 'f', s™^^) in the inclusive 
electron sample using the fitted branching fractions. These revised branching fractions are used to 
determine the b — > ceV background in the inclusive electron sample. 

The signal region > 2 GeV and s™*^^ < 3.5 GeV^/c^ was not examined until all selection 
criteria were fixed based on studies of simulated data and of control samples. Two control samples 
are used to reduce the sensitivity of the efficiency and background estimates to details of the 
simulation: (a) the DeV control sample described above; and (b) events satisfying the normal 
selection criteria but having s™'^^ > 4.25 GeV^/c'^ (high-s™^^ sideband), a sample dominated by 
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background. Efficiencies can be calculated, separately in data and MC, from the ratio of DeV 
candidates satisfying cuts (l)-(5) to the total number of selected DeV candidates having Eg, > 
2.0 GeV. The sensitivity of the h — > ueV signal efficiency to details of the simulation is reduced 
by multiplying it by the data/simulation ratio of DeU efficiencies. The high-s™'^'' sideband is used 
to normalize the Monte Carlo g™^ distribution to the data, reducing sensitivity to background 
normalization uncertainties. The total charmless semileptonic branching fraction is calculated as 



A' 



B{b 



rtiata 

.UK [ 



uev] 



^data 
9 ^MC De- 



N 



(2) 



BB 



where N^^^^ and -A''j^|j|*'* refer to the number of candidates in the signal and high-s™*^^ sideband 
regions of the data after subtraction of non-BB contributions determined on data taken below the 
BB threshold, and refer to background in the signal region and the yield in the sideband 

region in simulated events, e^^^ and are the efficiencies calculated on the corresponding DeV 
samples, pseudo-efficiency of the event selection determined from the data and Monte Carlo DeV 
control sample events, is the number of T{4S) BB decays analyzed. The total efficiency 

times acceptance for b — > ueV decays in the simulation is given by = €sig/« + ~ /«)' 

where /„ is the fraction of 6 — ueV decays generated in the signal region, and Cgig (e^) is the 
efficiency for an event inside (outside) the signal region to be reconstructed and pass our selection 
requirements. 

The effects of detector response and the boost of the B meson in the T(4S') frame are unfolded 
to produce a partial branching fraction that can be compared directly to theoretical calculations. 
We define 



AB{b^ueV) = B{b^ueV)U 



Ardata 
cand 



, .-^ Ardata 
a f MC ^side 

^"^bkg MMp- 
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data 
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jy-data _ A J-MC -^^idl^ 
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'BB 



(3) 



(esig/« + eiig(l -/«)) 



1 + 



ju 



- 1 



^sig 



-sig 



-1 



(4) 



Whereas the extraction of the total branching fraction B{b ueV) has a strong dependence on 
the signal modeling due to the theoretical uncertainty on /„, this dependence is suppressed in 
computing the unfolded partial branching fraction AH(6 ueV) since the ratio e^/egig ~ 0.029 
is much smaller than 1 (see Table 1). 

Figure 3 shows the electron energy and s^^^ distributions after cuts have been applied to 
all variables except the one being plotted. The yields and efficiencies are given in Table 1; these 
correspond to a branching fraction of ^2.37 it 0.22(stat)) x 10~^ using Eq. 2. We calculate the partial 
branching fraction AB for > 1.9 GcV in the B rest frame, which corresponds to E^ > 2.0 GeV 
in the T{4S) rest frame, and s™^^(/3 = 0) < 3.5 GeV^/c'*, unfolded for detector effects. We find 
ABib ^ ueV) = (4.51 ± 0.42(,tat)) x 10"^ 
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Figure 3: The electron energy distribution (top) and s^^^ distribution (bottom) after cuts have 
been apphed to all variables except the one being plotted. The arrows denote the signal region and 
also the high-s™*^^ sideband region (above 4.25 GeV^/c^). The number of background events from 
cascade decays and mis- identified electrons is small (events denoted as other). 



Table 1: Yields and efficiencies. All uncertainties are statistical except for on N^n. 



Data 


-^cand 


^side 


«(%) 




Nbb (10') 






8417 ± 164 


17776 ± 195 


8.72 ±0.52 




88.35 ± 0.97 


MC 


Mbkg 






esig(%) 




fu 




5687 ± 47 


17904 ± 81 


9.18 ±0.22 


3.26 ±0.03 


0.095 ±0.003 


0.1907 ±0.0002 



4 SYSTEMATIC STUDIES 

Systematic uncertainties are assigned for the modeling of the signal b — > ueV decays and background 
decays, and the modeling of detector response. The leading sources of uncertainty are listed in 
Table 2. Detector modeling for charged particle tracking, neutral reconstruction and charged 
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particle identification was evaluated by comparing data and MC in control samples. The simulated 
production rate of was verified by comparing data and MC momentum distributions for Kg in 
semileptonic decays, and the energy deposition of in the calorimeter was varied by a generous 
amount (±50%) to assess the corresponding uncertainty. The uncertainty due to Bremsstrahlung 
in the detector was based on the method used in Ref. [18]. QED final state radiation was simulated 
using PHOTOS [19]; comparisons with the analytical result of Ref. [20] were used to assess the 
systematic uncertainty. The uncertainty in the background was evaluated by varying the form 
factors and branching fractions of the b — > ceV decays. An additional uncertainty of 12% was added 
to account for the variation in the extracted branching fraction when the cut on was varied from 
1.8 GeV to 2.2 GeV. The modeling of signal decays is sensitive to the resonant structure at low mass. 
The exclusive branching fractions B{B — > hei'), where h = Tr,r), p,LO,ri' , were varied coherently by 
±30% to evaluate the uncertainty which proved to be negligeable. The sensitivity of the branching 
fraction to the parameters mf^ and a of the Dc Fazio-Ncubcrt model [13] was evaluated using 
uncertainties based on fits [14] to the b —>■ spectrum [6], leading to = 0.190 1q o24- 



Table 2: Systematic uncertainties on B{b — > ueu) and AB{b uev). 



Source of Systematics 


6B (%) 


5{AB) (%) 


la) Tracking efficiency 


±2.7 


±2.6 


2a) Electron ID efficiency 


+3.1 

-2.7 


+3.1 
-2.7 


3a) Charged particle ID 


+2.0 
-2.7 


+2.0 
-2.6 


4a) Bremsstrahlung 


+2.2 
-1.3 


+2.6 
-0.7 


5a) Neutrals reconstruction 


±6.3 


±6.4 


6a) Energy from 


+7.1 
-7.6 


+7.1 
-7.6 


7a) B counting 


±1.1 


±1.1 


A) Experimental systematics 


±10.9 


±11.0 


lb) B — ^ Xciv simulation 


+6.5 
-7.5 


±7.5 


2b) Radiative corrections 


+4.0 
-3.1 


+4.5 
-2.6 


3b) Stability scans 


±12 


±12 


B) Background simulation 


±14.4 


±14.6 


C) Signal simulation 


+25.4 
-12.5 


±4.6 


Total (A © B © C) 


+31.2 
-22.0 


±18.8 



5 RESULTS 

Using 88 million T(45') decays collected with the BABAR detector at the PEP-II e+e" storage 
ring, the charmless semileptonic branching fraction is determined using the method of Ref. [9]. 
The unfolded partial branching fraction is determined in the B rest frame for > 1.9 GeV and 
sgiax^^ = 0) < 3.5 GeV^/c^: 

AB{b ueV) = (4.51 ± 0.42 ± 0.50 ± 0.66 ± 0.19) x 10"^ (preliminary). 
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The uncertainties are from statistics (data and MC), detector modeling, background modeling and 
the modeling of 6 — > ueV decays, respectively. This partial branching fraction can be directly 
compared with theoretical calculations. We also determine 

B{b ueV) = (2.37 ± 0.22 ± 0.26 ± 0.34+° x lO'^ (preliminary), 

where the uncertainties are from statistics (data and MC), detector modeling, background modeling 
and the SF parameters and a, respectively. The uncertainty due to the modeling of resonant 
states in 6 ^ ueV decays is negligible compared to the other uncertainties. 
We extract from our measurement using [21] 



iKbl = 0-00424^ ^2q2 ' tb ^ ^'^'^^^^^ ^ 0-039i/^3), 

where the first coefficient and the uncertainties quoted in Ref. [21] have been updated using the ex- 
perimental input from the moment measurements obtained in Refs. [18, 22, 23]. The first coefficient 
takes also into account electroweak radiative corrections [24]. 
Taking tb = 1.604 ± 0.012 ps from [25] we find 

= (4.57 ± 0.21 ± 0.25 ± 0.34 t^fg ± 0.22) x 10"^ (preliminary), 

where the first four uncertainties are as for the total branching fraction and the last comes from 
the HQE relating \Vub\ to the full branching fraction. 

The theoretical uncertainty is dominated by the SF parameter . The variation taken in this 
analysis is large [14]: = 4.735loi55 GeV/c^. A more precise estimate of the SF parameters has 
been determined from the B Xg^j photon energy spectrum measured by Belle only recently [26] 
from which we obtain the following preliminary results: 

AB(6 ^ ueV) = (4.46 ± 0.42 ± 0.49 ± 0.64 ± 0.19) x 10"^, 

B{h ueV) = (2.76 ± 0.26 ± 0.30 ± 0.40;[|:^||;[|:^^) X 10"^ (/« = 0.1628 ± 0.0002), 

where the uncertainties arc from statistics (data and MC), detector modeling, background modeling, 
the SF parameters m^^ and a, and the uncertainty due to the modeling of resonant states in h ueU 
decays, respectively, and 

iKbl = (4.99 ± 0.23 ± 0.25 ± 0.34 1°;°^ ± 0.22) x 10-^ 

where the first five uncertainties arc as for the total branching fraction and the last comes from the 
HQE relating \ Vub\ to the full branching fraction. Note that all of the above results are preliminary. 

Current theoretical work (e.g., see Ref. [27]) is aimed at connecting the parameter in the 
SF with determinations of from moments in 6 — ceV decays or from the ^(15) mass, and 
may significantly reduce this uncertainty. The results obtained in this analysis have very little 
dependence on the SF parameter a and are complementary to those obtained from studies of the 
electron endpoint and of the mass of the recoiling hadron in semileptonic decays. 
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